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J. MyslivecekQ and T. Jarolfmek 

Department of Electronics and Vacuum Physics, Faculty of Mathematics and Physics, Charles University, V Holesovickdch 2, 

180 00 Praha 8, Czech Republic 



<3\ 



> 

o 
in 
o 

On 

> 



c 

o 
o 



X 



P. Smilauer 

Institute of Physics, Academy of Sciences of the Czech Republic, Cukrovarnickd 10, 162 53 Praha 6, Czech Republic 

B. Voigtlander and M. Kastner 
Institut fur Grenzflachenforschung und Vakuumphysik, Forschungszentrum Julich, 52425 Jiilich, Germany 

(February 1, 2008) 

Surface reconstructions can drastically modify growth kinetics during initial stages of epitaxial 
growth as well as during the process of surface equilibration after termination of growth. We 
investigate the effect of activation barriers hindering attachment of material to existing islands 
on the density and size distribution of islands in a model of homoepitaxial growth on Si(lll)7x7 
reconstructed surface. An unusual distribution of island sizes peaked around "magic" sizes and 
a steep dependence of the island density on the growth rate are observed. "Magic" islands (of 
a different shape as compared to those obtained during growth) are observed also during surface 
equilibration. 
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Investigation of island structures formed during the 
initial stages of epitaxial growth allows us to explore ki- 
netic mechanisms that govern the ordering of deposited 
atoms.El A lot of attention has recently focused on the 
time- and growth-conditions dependence of the island 
densitya as well as on the bell-shaped distribution of is- 
land sizesa whose origin can be traced back to the dis- 
tribution of island capture zones. However, real growth 
systems are invariably more complicated than the ide- 
alized models of epitaxy commonly used. For example, 
the presence of surface reconstructions can completely 
change the growth behavior. 

In homoepitaxy of Si on Si(lll)7x7 reconstructed sur- 
face, a process of "reconstruction destruction" was de- 
scribed by Tochihara and Shimada.Q The need to can- 
cel surface reconstruction around a growing island gives 
rise to barriers to attachment of new material to existing 
islands. Growth with barriers to attachment has been 
already studied theoretically: The dependence of the is- 
land density on srawth conditions was explored using 
analytic methods]^ while kinetic Monte Carlo (KMC) 
simulations of a simple growth model revealed an island 
size distribution multiple-peaked around "magic" sizes .0 

Here we present a detailed KMC model of 
Si/Si(lll)7x7 molecular beam epitaxy (MBE), with 
barriers to attachment included. With the help of this 
model we investigate the time- and growth-rate depen- 
dence of the island density, the shape of the island-size 
distribution, as well as island decay and filling of vacancy 
islands on the surface. The results of our simulations 
compare favorably to available experimental data about 
the Si/Si (111)7x7 system. We also discuss those fea- 
tures of the model kinetics that are specific to growth 
with barriers to attachment. 

Dynamics of Si/Si(lll)7x7 MBE growth was experi- 



mentally studied by Voigtlander et aZ.i~0 The most in- 
teresting feature observed, the existence of kinetically 
stabilized magic sizes in the island size distribution, was 
reported and KMC modeled in Ref. |7[ The model dis- 
cussed in this work is a generalization of the model from 
Ref. |t] which is in turn based on the "reconstruction de- 
struction" model of Si/Si/lll)7x7 growth proposed by 
Tochihara and Shimada.u The material is deposited in 
units (cf. below) that are randomly placed onto sites of 
a honeycomb lattice. These sites represent half-unit cells 
(HUCs) of Si(lll)7x 7 surface reconstruction. Two types 
of HUCs exist: In unfaulted HUCs (marked U), a sur- 
face atom bilayer follows bulk bilayer stacking, in faulted 
HUCs (marked F), the surface bilayer is 30° rotated with 
respect to the bulk, forming a stacking fault O Material 
in a HUC may be either non-transformed (marked X, 
models the free Si adatoms diffusing on the reconstructed 
Si surface) or transformed (marked T, corresponds to a 
HUC, where the material underwent the reconstruction 
destruction process and subsequent crystallization). The 
model thus becomes de facto a two-species one. The non- 
transformed units of the material randomly walk on the 
surface, meet each other and transform. Transformation 
is an activated process, activation energy for transforma- 
tion is supposed to be higherJn F HUCs due to the need 
to remove the stacking fault .□ 

In the experiment, deposited material does not diffuse 
as HUCs£3 The use of units of deposited material al- 
lows us to model easily collective processes during "re- 
construction destruction" around island edges. Since, in 
general, the processes at island edges determine behav- 
ior of growth models (as they determine growth behavior 
of real systems), the simplification of material deposi- 
tion and surface diffusion should not affect the model 
behavior in a substantial way. Our simulation scheme 
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is coarse-grained and ignores all processes on the length 
scales smaller than a HUC and on the time scales shorter 
than the time required to transport material from one 
HUC to its nearest neighbor. 

The hopping rate for a HUC hop is 
VD=voCx.p(-E D /k B T) where E D =E S + (x + t)E§ for 
an X-HUC, E D =E S + xE§ + tE% for a T-HUC, x and 
t being the numbers of X and T neighbors, respectively, 
E^ bond strengths of X — X and X — T pairs, Ejj 
bond strength of a T — T pair, and Eg the surface bar- 
rier to diffusion. The rate of an X-HUC transformation 
is vt=vq exp(— Ex/hsT) where Ej, = Ea — tE e d ge for 
an F-HUC, E% = E A - tE cdgc - E dm for a U-HUC, 
Ea being the barrier to attachment, E cdgc a decrease 
in the barrier due to a transformed neighbor, and E di g 
the barrier difference for F- and U-HUC overgrowth. 
Transformation of an island begins at an F-HUC with 
>2 X-neighbors, and the rate of this nucleation process 
is v=vq exp[(£r — E cdge )/kBT]. The model has seven 
parameters, here we report results for fo=10 13 s" 1 , 
E s =1.5 eV, #£=0.3 eV, £$=0.1 eV, E A =2.3 eV, 
-E'edge=-E ; diff=0.35 eV, which gave the best agreement 
with experimental results. Using the model, we tried to 
reproduce both growth and equilibration processes on the 
Si/Si(lll)7x7 surface on a real time- and spatial scale. 
The HUC in the model is thus considered to be 1 bilayer 
(BL) of Si(lll) thick, of a triangular shape with the 
edge length of a— 26.9 A and consisting of 49 Si atoms. 
All calculations were performed on 200x200 HUC lattice 
with periodic boundary conditions. 

Growth. We tried to fit the exponent x given by 

N^F X (1) 

where the dependence of the island densi£ji,iV on flux F 
at a constant temperature is measured. El'EiEl The experi- 
mental value of x=0.75 for T=680 K and T=770 K was 
reported in Ref. |l3"[ In this experiment, samples were pre- 
pared at a given temperature by deposition of «0.15 BL 
Si on the Si(lll)7x7 surface followed by rapid quenching 
to room temperature. The experimental morphologies of 
the layer can therefore be considered snapshots from the 
Si/Si(lll)7x7 surface morphology evolution. 

In the model, we calculated the flux dependence of the 
density of transformed (i.e., crystalline) islands N at con- 
stant total coverage Otot (Ref. [l4]). Results are shown in 
Fig. [I]. At lower fluxes, a power-law N=F X dependence 
with X6so=0.76±0.03, X77o=0.75±0.04 is observed. In 
the experiment, disordered growth occured at high fluxes. 
In the model, deviations from the power-law behavior of 
N=N(F) are observed. 

The presence of the barriers to attachment in the 
model introduces specific features into its dynamical be- 
havior. In the dependence of N on the total coverage 
©tot (Fig. |^a) the nucleation onset (fast increase of N) 
takes place at higher coverages than in standard (mini- 
mal) growth model (Fig. |b), Ref. [if. The maxima of 
A(Otot) are sharp and occur at higher to t than the 
broad maxima of the standard model. 



Island growth in our model starts by a transformation 
of 3 neighboring HUCsD Due to the barriers to attach- 
ment this happens only after a certain time elapsed from 
HUC clustering . This time does not scale with flux 
in contrast to the time for clustering of units of mate- 
rial. Transformation becomes the rate-limiting process 
of island formation. By the time the transformation 
begins, more material has been deposited for the case 
of higher flux F. This shifts the nucleation onset and 
A max (8 to t) position to higher 9 to t- With increasing F, 
amount of non-transformed material at a given Otot in- 
creases (Fig. [2)3) , and so does the disorder of the quenched 
sample. 

A newborn island composed of 3 units can further 
grow or decay. The competition between growth and 
decay of transformed clusters results in a steep decrease 
of N(Q tot) for N > A max . This decrease is not due to 
coalescence, because the islands at A=A max are small. 
Both the delay in island formation and the instability 
of islands leave traces in the time evolution of the mean 
island size, cf. Fig. ||[ 

The shift of A max to higher O to t than in standard 
model and the fast decrease of A(0 to t) after reach- 
ing A max contribute to the high value of \ observed 
at Gtot=const during growth with barriers to attach- 
ment. With the barrier to attachment decreasing and the 
nearest-neighbor bond strength increasing, x decreases. 

Analytical theories usually relate the value of the 
growth exponent ^ to i*, the number of material units 
in a "critical" (i.e., largest unstable) island in the growth 
system. For the determination of i* , we cap use a for- 
mula x = 2i*/(i* + 3) derived by Kandel.EE3 x m tne 
model varies smoothly with E$, E*, and Ejj so that 
the corresponding values of i* are non-integer numbers 
between i*=l (x=0.5) and i*=2 (x=0.8). In addition, 
a dependence of the island density on E$ (the substrate 
contribution to the hopping barrier) was found, in con- 
trast to predictions of Ref. g. 

The instability of growing islands does not strongly 
affect morphologies obtained for Si/Si(lll)7x7 MBE 
growth. We observe the characteristic multiple-peaked 
island size distribution (Fig. ||c, d), but with broader 
peaks (in better agreement with the experimental results) 
as compared to the model with detachment of material 
from islands forbidden.Q In the morphologies of experi- 
mental samples, non-transformed clusters of Si adatoms 
formed during quenching of samples are visible (Fig. |]b). 
The density of non-transformed material may be experi- 
mentally determined and compared to the model results. 

Equilibration. The decay of Si/Si(lll)7x7 recon- 
structed 1 BL-high adatom (A) and vacancy (V) islands 
was experimentally studied in Refs. |l7j-|l9|. Using STM 
at an elevated temperature, the authors of Refs. [l7| and 
|l8| followed a number of isolated (a nearest island or a 
step edge at a distance more than 800 A) adatom or va- 
cancy islands and studied the temperature dependence 
of their decay rates. The decay rates v A , v v showed 
Arrhenius behavior v a > v = V q' v exp(-E A ' v /k B T) 
with ^=2xl0 11±1 adatoms-s- 1 , £^=1.5±0.1 eV for 
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adatoms, 



=3xl0 9±1 adatoms-s-\ £ y =1.3±0.2 eV 



for vacancies, respectively. The decay rate of vacancy 
islands was found to be approx. 5 times lower than that 
of adatom islands .E3 

With our model, we traced the evolution of a 96-HUC 
compact adatom or vacancy island placed on a vicinal Si 
surface (U-type steps, the terrace width of 480 nm, the 
distance from the descending step edge of 240 nm) equili- 
brated at a given temperature. Step edges on the vicinal 
surface form adatom sources and traps necessary for true 
disappearance of a single adatom or vacancy island in a 
model with periodic boundary conditions. 

The temperature dependence of the decay rates for 
adatom and vacancy islands in our model is shown in 
Fig. [5^,, b. With the parameters listed above, the de- 
cay rates in the model are higher than the experimen- 
tal ones (i/^=3xl0 14 - 5±0 - 5 adatoms-s~ 1 , £^=2.1±0.1 eV, 
^=lxl0 14±1 adatoms-s-\ £^=2.0±0.1 eV), and the 
decay rate of vacancy islands is approx. 2 times lower 
than the decay rate of adatom islands. 

The authors of Ref. 17 18 attributed the difference be- 
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tween the decay rates of adatom and vacancy islands to 
the effect of the Ehrlich-Schwoebel (step-edge) barrier in 
the Si/Si(lll)7x7 system. We do not believe that the 
Ehrlich-Schwoebel barrier plays any role: Growth exper- 
iments provide no compelling evidence of the presence of 
an appreciable Ehrlich-Schwoebel barrier at step edges 
on Si/Si(lLl,)7x7 surface within the relevant tempera- 
ture range. £3 

We also modeled adatom and vacancy islands 
decani with the barrier to attachment "switched 
ofp'Jij The decay rates thus obtained were 
lower (^=5xl0 13±1 adatoms-s-\ £^=2.1±0.2 eV, 
^=lxl0 12±1 adatoms-s- 1 , £^=2.1±0.2 eV), but the 
decay rate of vacancy islands was still approx. 2 times 
lower than for adatom islands. This observation agrees 
with results of a standard growth model on square 
lattice.Eil The difference in these decay rates on the vici- 
nal surface thus seems to originate from the difference, of 
geometry of adatom and vacancy island boundaries. Ell 

In Fig. [|c, a typical time evolution of the size of a de- 
caying island in a model with the barriers to attachment 
is shown. We see that stable ("magic") Si islands do 
exist. They correspond, to, equilibrium island shapes ex- 
perimentally observedt3i!3 and differ from magic shapes 
observed during Si/Si(lll)7x7 growthE Magic islands 
are compact (2 nearest-neighbors for all perimeter HUCs) 
and the barrier to attachment prevents their shape from 
"being spoiled" by attachment of material surrounding 
the island. No stable shapes are observed during island 
decay for the model without barriers to attachment. 

In this work, we presented a coarse-grained model of 
Si/Si(lll)7x7 MBE growth with an activation barrier 
to attachment of new material to existing islands imple- 
mented. We demonstrated that this barrier contributes 
to the steep growth-rate dependence of the island den- 
sity observed in Si/Si(lll)7x7 MBE and helps to stabi- 
lize "magic" island shapes in both growth and relaxation 
experiments. 
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FIG. 1. Flux dependence of the island density N for 
Si/Si(lll)7x7 MBE is N=F X with x=0.75. Experimental (o) 
(Ref. |l|) and modeled (•) N(F) dependencies for 680 K (a) 
and 770 K (b) are shown. In the model, N at 6 t ot=0.15 BL 
was measured. 
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FIG. 2. Evolution of the island density vs. total coverage 
for a model with the barriers to attachment depends on flux 
(a) and differs from that of a standard growth model (b), 
Ref. [l5|. The barrier to attachment limits the rate of creation 
of reconstructed islands. At high fluxes, most of the material 
present at the surface is non-transformed (c). Data for (1) 
F=10- 4 BL-s-\ (2) F=Kr 3 BL-s-\ and (3) F^lO^BL-s" 1 
are shown. 
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FIG. 3. The delay of island creation is clearly visible in the 
record of the mean island size time evolution. Instability of 
newborn islands causes the decrease of (s) at short times. 
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FIG. 4. Growth morphologies in the model (a) and exper- 
iment (b). During growth, magic island sizes are stabilized, 
resulting in a non-trivial island size distribution (c, d), Ref. ^. 
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FIG. 5. Arrhenius plots of experimental (o) (Ref. |l7| ) and 
modeled (•) decay rates for adatom (a) and vacancy (b) is- 
lands. In the model, the rate of vacancy filling is 2x lower 
than that of adatom islands decay. During island decay, 
"magic" island sizes are stabilized (c). Insets show some of 
the observed stable island- morphologies. These are close to 
equilibrium island shapesE-jlij and differ from "magic" island 
shapes observed during growth (cf. Fig 0b) . 
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